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ABSTRACT 

In  this  report  we  discuss  adaptive  array  per  foresee  when  a 
continuous  sector  of  interference  is  incident  upon  the  array.  A 
model  is  developed  for  4 two-element  array  with  «n  arbitrary  number 
of  incident  signals.  From  this  mode)  the  array  weights  are  derived 
and  formulas  for  the  array  performance  are  developed.  It  is  then 
shown  how  this  model  is  readily  extended  to  the  case  of  a continuous 
interference  sector. 

from  these  analytical  results,  the  array  performance  is  then 
examined  for  specific  cases  of  spatially  dispersed  interference.  It 
is  seen  that  the  array  performance  is  not  always  degraded  as  the 
interference  angular  sector  becomes  larger. 
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I.  INTRODUCTION 

Signal  processing  antenna  arrays  have  been  studied  for  many 
years.  For  receiving  arrays  used  in  communication  and  radar  systems, 
it  Is  important  that  the  output  signal  to  noise  ratio  (SNR)  be 
optimized.  It  is  highly  desirable  that  this  optimization  be  auto* 
matlc.  That  is.  for  each  change  in  the  antenna  or  signal  environment 
(e.g.,  change  in  direction  of  arrival  of  the  signals,  change  in 
electronic  components,  etc.),  the  system  should  modify  itself  until 
the  SNR  is  again  optimum.  One  method  of  achieving  this  modification 
is  through  internal  feedback.  Antenna  systems  of  this  type  are  called 
adaptive  arrays. 

Shor  [Y]  was  one  of  the  first  to  suggest  an  adaptive  process 
which  maximized  the  SNR  of  an  array  of  hydraphones.  Wldrow  fi]  et 
aj_  suggested  a feedback  technique  which  minimized  the  me  an -square 
error  between  the  array  output  and  a reference  signal.  Applebaum  |3j 
and  Griffiths  f4)  discussed  similar  concepts.  An  early  experimental 
adaptive  array  was  built  by  Riegler  and  Conpton  [5] . 

The  behavior  of  an  adaptive  array  with  spread  spectrum 
communication  signals  has  been  studied  by  Relnhard,  Huff,  Conpton 
and  others  [7,8,10,11,13,14} . Beml  [9]  has  suggested  a method  of 
angle  of  arrival  estimation  using  an  adaptive  array.  A four-element 
array  capable  of  arrival  angle  estimation  and  sensor  comnunl cations 
has  been  implemented  and  experimental  results  are  described  by 
Swamer  and  Berni  [15]. 


In  this  report  we  study  the  performance  of  adaptive  arrays 
based  on  the  IMS  algorithm  [2,3,4]  when  interference  is  incident  on 
the  array  from  an  angular  sector  of  finite  width.  Previous  studies 
have  concentrated  on  the  case  of  interference  signals  that  arrive  from 
a single  direction  is  space,  he  begin  by  assuming  a large  number  of 
interference  signals  to  be  incident  on  the  array  (many  more  than 
degrees  of  freedom  of  the  array).  The  signals  are  assumed  to  be  un- 
correlated and  to  arrive  from  different  direction*  within  a certain 
sector.  The  array  performance  is  studied  as  a function  of  the  size  of 
the  sector.  The  case  of  discrete  Interfering  signals  is  generalised 
to  the  case  of  radiation  arriving  from  a continuous  sector. 

In  section  II  certain  mathematical  preliminaries  are  investi- 
gated. By  analyzing  a one-element  array,  the  effect  of  time  vary- 
ing coefficients  in  the  differential  equations  for  the  weights  is 
studied.  It  is  found  that  under  certain  conditions  these  time  varying 
components  may  be  neglected  without  greatly  affecting  the  weight 
solutions.  In  section  Itl  the  general  solution  for  the  weights  of  a 
two-element  array  with  an  arbitrary  number  of  incident  CW  signals  is 
developed.  In  sections  IV  and  V this  solution  is  used  to  study  the 
system  performance.  In  section  VI  the  case  of  many  discrete  signals 
is  extended  to  a continuous  column  of  inpinging  radiation.  The 
weight  solutions  are  derived.  Finally,  in  section  VII  numerical 
examples  of  array  performance  in  a many-signal  environment  are  presented. 
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11.  THE  FEE06ACK  ALGORITHM  AMO  THE  ARRAY  STRUCTURE 

An  N-element  antenna  array  with  quadrature  weighting  of  the 
Input  signals  Is  shown  in  Fig.  1.  In  this  figure  the  quadrature 
signals  from  the  first  antenna  are  denoted  *j(t)  and  *,(t),  those 
from  the  second  are  *,(t)  and  xfc(t),  and  so  forth.  The  weights  asso- 
ciated with  these  2N  signals  (N  antenna  elewmts)  are  similarly 
Indexed.  The  array  output  may  be  written  as 

2H 

(1)  *(t)  • l WjXjft)  . 

1-1 

We  define  yj(t)  to  be  the  signal  incident  an  the  1-th  antenna.  xf(t) 
will  be  referred  to  as  the  eleeent  signal.  The  error  signal  is 

2 N 

(2)  c(t)  » R(t)  - l w^ft) 

1»1 

where  R(t)  is  the  reference  signal.  Realistically,  it  must  be  assumpd 
that  the  desired  signal  contains  modulation  components  that  are  u»- 
known  at  the  receiver.  Hence,  the  reference  signal  cannot  be  made 
exactly  equal  to  the  desired  part  of  the  Incoming  si yral,  but  can  only 
approximate  It  In  some  sense.  For  correct  operation  In  an  inter- 
ference rejection  system.  It  Is  necessary  to  generate  a reference 
signal  which  resembles  the  desired  signal  and  correlates  poorly  with 
interference  [1425  |.  The  array  will  then  act  to  drive  the  error 
signal  to  zero  in  the  mean-square  sense.  Such  Interference  rejection 
is  achieved  by  placing  spatial  nulls  In  the  direction  of  interfering 
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transmissions.  One  way  to  obtain  a reference  signal  is  to  derive  it 
from  the  array  output.  Certain  difficulties  may  arise,  however,  in 
attempting  to  derive  the  reference  sicyial  from  the  array  output.  Some 
of  these  are  discussed  In  detail  In  [8,14,25]. 

From  Eq.  (2)  the  squared  error  becomes 

2N  2N  2N 

(3)  c*(t)  - a7(t)  - 2R(t)  l w x ft)  ♦ l l w.w.x.(t)  x (t). 

1-1  11  1-1  j-1  1 J 1 3 

The  mean-square  error  is  thus: 

2M  2N  2N  

(4)  c*(t)  - R*(t)  - 2 [ wj  R(t)x,(t)  ♦ l l w1wixi(t)x1(t) 

1-1  1-1  j«l  J J 

where  the  overbar  represents  the  action  of  a low-pass  filter  as  will 
be  discussed  later  In  this  section. 

Differentiating  Eq.  (4)  with  respect  toWj  yields 

(5)  * [JitT]  • - -2x1  (t)e(t)  . 

1 ® J 

The  feedback  system  Is  based  on  the  so  called  IMS  algorithm 
[2,3,4  ].  Each  veight  Is  controlled  by  the  relation 

<«>  5T  ' ' [t2(t)]  • 

Then  from  Eq.  (5),  the  feedback  equation  becomes 
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dw. 

(7)  ~ - 2k  c(t )xAt) 

dt  ' 

which  leads  to  the  feedback  structure  shown  in  Figs.  2a  and  2b.  Wh^n 
Eq.  (2)  is  used  to  substitute  for  t(t)  In  Eq.  (7),  and  all  the  terms 
involving  are  collected  on  the  left,  it  is  fouid  that  the  weights 
satisfy  the  system  of  differential  equations  given  by 

dw.  2N  ______ 

(8)  — 1 ♦ 2k  [ [xt(t)xt(t)]w.  • 2k[R(t)xi(t)]  . 

dt  jBl  J J 

The  system  of  Eq.  (8)  may  be  solved  (under  certain  simplifying  ( 
conditions)  to  yield  the  time  response  of  the  weighting  coefficients 
Wj(t).  From  the  weight  solutions,  the  performance  of  the  system 
under  various  conditions  may  be  studied. 

Before  examining  the  method  of  solution  we  discuss  the  meaning 
of  the  overbar.  We  have  stated  previously  that  the  overbar  represents 
the  action  of  an  ideal  low-pass  filter.  For  example.  If  x^(t)  and 
R(t)  in  Eq.  (4)  each  contain  two  CW  signals  of  different  frequency, 
their  product  will  contain  d.c.  terms  along  with  components  at  the 
sum  and  difference  frequencies.  Typically,  multipliers  used  in  an 
adaptive  array  will  not  pass  the  sum  frequencies.  For  example,  in 
one  array  implemented  at  Ohio  State  [15],  the  array  processing  was 
done  at  70  M-Hz.  Trans -conductance  multipliers  were  used,  which  did 
not  pass  frequencies  higher  than  100  M-Hz.  Therefore,  sum  frequency 
terms  (at  140  M-Hz)  were  not  passed. 

6 


The  difference  frequency  terms,  however,  were  within  the  pass- 
band.  In  previous  studies  of  adaptive  arrays  with  CW  signals,  the  effects 
of  these  difference  frequency  terms  have  been  neglected[5-ll ,1 3,15,16 ,21 -25]. 
In  this  section,  we  discuss  the  conditions  under  which  these  terms  moy  be 
neglected  without  affecting  the  solutions  of  the  weights.  It  will  be 
shown  that  there  exists  a cutoff  frequency  beyond  which  the  effect  of 
these  difference  terms  may  be  neglected.  It  will  also  be  demonstrated 
that  the  power  contained  in  frequency  components  below  this  cutoff  fre- 
quency will  be  small  compared  to  the  power  in  the  d.c.  component.  These 
conditions  will  allow  us  to  neglect  all  but  the  d.c.  terms  in  Eq.  (8). 

Let  us  suppose  that  a one-element  array  has  two  CW  signals 
incident  upon  it.  The  signal  is 

(9)  y(t)  =V2a[  cos^t)  +JI  a p cos(u»2t). 

We  assume  ideal  reference  signal  processing  is  available 
(i.e.,  the  reference  signal  is  perfectly  coherent  with  the  desired 
signals  and  is  unaffected  by  the  weights),  and  the  reference  signal 
R(t)  is  given  by 

(10)  R(t)  = cosUjt)  + cos(w2t). 

In  other  words,  R(t)  is  coherent  with  both  signals  (i.e.,  we  have 
chosen  both  signals  to  be  desired).  Upon  quadrature  splitting  we 
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obtain  for  the  element  signals 

(11)  x (t)  • a cos (u  t)  ♦ a cos(w  t) 

* * l 2 2 

and 

(12)  x2(t)  s sinfw^)  + a si n(w2t). 

Since  the  cverbar  represents  an  ideal  low-pass  filter  with  a 
cutoff  frequency  higher  than  Ul-u,2  but  lower  than  2^  and  Zu  , the 
terms  Xjlt)Xj(t)  are 

ai  + 

(13)  (t)Xj (t)  = - + flla2  cos[(Ui-w2)t] 

a2  + a2 

(14)  x2(t)x2(t)  = — 2 2 + a^  cos[(Ui-<o2)t] 

and 


(15)  x](t)x2(t)  = x2(t)Xi(t)  = 0. 


1 

i 

> 


We  also  find 


(16) 


R(t)Xj (t)  = 


31  + fl2 


ai  + d2 


C0S[(Wj-a»2)t], 
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and 


(17) 


RTtTxJtT 


♦ a- 


sin[(u,-w.  )t] 


When  Eqs.  (13)  - (17)  are  substituted  Into  Eq.  (8),  the 
differential  equations  describing  the  weights  can  be  written  in 
matrix  form  as 


/w,(t)> 

( 

t\w2(t) J 

) + 2k 

V 

»'l+a2 

+ a^cosfAwt) 

V 0 


\ 


a‘?+a‘ 


) a2 

^ + ala?  COS(Awt) 


|C/,(t) 

w?(t) 


( a !+a2  + (a 1+a2)cos(Aut)^ 
-(a1+a2)sin(Au)t) 


where 

(19)  Aw  = W]-u? 

Let  us  examine  the  equation  for  weight  Wj(t).  We  notice  that  the 
system  of  Eqs.  (18)  is  uncoupled,  so  we  can  solve  for  the  weights 
independently  of  one  another.  The  equation  describing  the  reponse  of 
the  in-phase  weight  (Wj(t))  is  then 

dw i ( t ) 

(20)  — rz + [A  + B*cos(Au>t)]wi(t)  = [C  + D*cos(Awt)]  u(t) 
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where 

(21)  A-  k •(•}♦•!) 

(22)  B = 2 • k • a,  • a , 

(23)  C = k • (a,  + a2) 

(24)  D = k • (aj  + a2) 

and  u(t)  is  the  unit  step  function.  We  arbitrarily  choose  the  initial 

value  ofwx(t)  to  be  zero  at  t = 0-.  The  Fourier  transform  ofw,(t) 
exists  and  Eq.  (20)  may  be  written  in  the  frequency  domain  as 

(25)  (jw)  • Wjfw)  + A * Wjfu)  + j • [Wj  (w+Au>)  + Wj  (w-Aoj)  ] = 

C * [tt  * 6(u>)  + i—  ] + • [6(u-A<o)  + 6(u>+A«)]  + ^ 

J“  2 (A«)2-«2 
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where  W[(w)  is  the  Fourier  transform  of  wj(t),  6(w)  is  the  Dirac 
delta  function  and 

(26)  j - p\. 

Solving  Eq.  (25)  for  Wj(w)  and  taking  the  inverse  Fourier 
transform  shows 

(27)  W,(t)  = [e"At  • u(t)3*  [C  • u(t)]  + [e'At  • u(t)]  * [D  • cos(^t) 


[B  • Wj(t)  • cos(Awt)]  * [e'At  • u(t)] 


where  represents  convolution  in  the  time  domain,  defined  by 


(28)  f(t)  * g(t)  = f f(t-x)  • g(x)dx  = f f(T)  • g(t-x)dx  = F(w) 


The  first  two  convolutions  of  Eq.  (27)  are  easily  performed.  The 
third  may  tie  written  in  integral  form  yielding  as  an  expression  for 

Wj(t): 


(29)  Wj(t) 


. C C -At 
" A • A e 


+ 


A*  D>e 


-At 


A ?-+( Aw)2 


i D* rA«cos(flwt)-*-(Ah)) -sin(Au)t)  1 
A2  + (M2 


• u(t)3  I 


G(“  ) • I 


j 

I 

| 
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COS(Awt)  • w,(r)di  . 


B 


rt  -”A(t-T) 

I e • 

* -'H 


We  can  generate  an  asymptotic  expansion  for  *j(t)  valid  for 
large  (au).  Continually  Integrating  the  final  term  of  Eq.  (29)  by 
parts  will  yield  a solution  of  the  form 


(30) 


B 


e A^t‘T^cos(/w»T)*w,(T)*dT 


J w,(t)  ♦ 

A2  ♦ (aw)2 


a,(t,A«..) 

w.(t) 

(A2+(a*)2)2  * 


a 3(t,Aw) 

(A2+( Aw)2) 3 


»,(t) 


When  this  result  Is  used  In  Eq.  (29),  the  solution  for  w,(t)  may  be 
approximated  by  the  first  term  of  the  series  if  Aw  Is  high  enough. 
Integrating  Eq.  (30)  by  parts  we  find 


(31)  B 


e T^*cos(AuT)*Wj(T)dT 


£A-B«cos(Ato)t)+B-  (Awl-sinlAutll 

A2  ♦ (Aw)2 


Wj(t) 


B*e 


■At 


A"+( A u)‘ 


eAx(A 


•COS(Awi)  + (Au))*sin(A(DT))w1(T)dT  . 


Recall  from  Eq.  (20)  that 
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(32)  Wj(t)  * [C+D-cos(Awt)]u(t)  - [A+B’COs(Awt)]  • wj(i). 

Substituting  for  w(T)  Into  the  last  Integral  of  Eq.  (31)  only  those 
terms  of  Eq.  (32)  which  contribute  to  aj(t,Au>)  (1,e.»  wt(T)  ■ 

D cos(au>t)  u(t)  - B cos(awt)  Wj(t))  yields 

(33)  — - [ e^T(A‘cos(Awx  )+(a«) *sin(A«i»T ) ) 'w.  (t )dis 

A2+(Aw)2  1 


B-e~At  f 
A2 + ( Au> ) 2 ■'0 


(A0-A*B'Wj(T))dT  . 


Continued  integration  by  parts  of  Eq.  (34)  again  results  in 
integrals  containing  *j(t).  Substituting  in  each  case  for  *,(t) 
only  those  terms  of  Eq.  (32)  which  contribute  to  a^t.Aui)  enables  us 
to  write  Eq.  (33)  as 


(34) 


B«e~At  ft 
A2+(au)2  *0  2 


(A.D-A*B’W2(T))dT 


B-D 

A2+(Aw)2 


(l-e"At) 


~B  + B •C»t-e~At 

A2+(a«)2  A2+  (a w)2 
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We  can  then  substitute  the  result  of  Eq.  (34)  for  the  integral 
expression  in  Eq.  (31).  Upon  substituting  Eq.  (31)  into  Eq.  (29),  we 


(35)  w (t)  * j-t  e'At  +~'D,eV+  - 

* M a*  , / . \Z 


A2-f(Aa))2  A?  ♦ (Aw)2 


(A*B*cos(Awt)+B- (Aw)*sin(Awt))*w  (t)  B»D«(l-e~At) 


A2  + (Aw)2 


A2  ♦ (Aw)2 


+ (Aw)2  A2  + (Aw) 


We  can  now  solve  Eq.  (35)  forw](t).  After  a long  period  of  time 
(i.e.,  after  the  transients  have  died  out),  we  obtain  forWj(t): 
(36) 

c A-D-cos(Awt)+D-(Aw)-sin(Awt)-C-B-cos(Awt)  - — 
w,(t)  ^ + — — A 

t->0o  A2+  (Aw)2+  A-B*cos(Awt)+B'  (Aw)*sin(Awt)  + B2 


The  first  term  on  the  right  side  of  Eq.  (36)  is  the  steady-state  term 
which  would  occur  if  B*cos(Awt)  and  D«cos(Awt)  were  neglected  in 
Eq.  (20).  Eq.  (36)  is  bounded  by 


15 


(37)  w (t)  < £ ♦ A‘D  * 

* A A ♦(Aw)  +B  -AB-8-(Aw) 

In  other  words  we  have  replaced  the  sin  (Awt)  and  cos(Awt)  terms  with 
unity  in  the  numerator  and  minus  one  in  the  denominator.  From  Eq.  (37) 
we  notice  that  the  second  term  will  be  small  compared  to  C/A  when 

(38)  Aw  >>  = k • (a^  ♦ ap  = wQ. 

In  other  words,  when  the  difference  frequency  Aw  is  much  larger  than 
the  product  of  the  loop  gain  constant  (k)  and  the  sum  of  the  signal 
powers,  the  steady-state  solution  for  Eq.  (20)  will  be  approximately 
C/A.  We  see  then,  that  when  the  inequality  given  in  Eq.  (38)  is 
satisfied,  the  difference  frequency  terms  (B*cos(Awt)  and  D*cos(Awt)) 
in  Eq.  (20)  may  be  neglected.  For  example,  in  one  adaptive  array 
built  at  Ohio  State  [15 j , the  value  of  w0  was  typically  between 
100  Hz  and  1 k-Hz.  In  this  array,  difference  frequencies  greater 
than  lk-Hz  did  not  affect  the  weight  solution. 

In  an  adaptive  array  for  conmuni cations,  the  signal  x^t)  will 
be  a bandlimited  process  centered  at  some  non-zero  carrier  frequency. 
(For  exanple,  the  Ohio  State  adaptive  array  [15 | operated  at  70  M-Hz 
with  a 10  M-Hz  bandwidth.)  It  will  be  adequate  for  the  present 
discussion  to  assume  that  x^(t)  has  a flat  power  spectral  density  of 
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height  a and  bandwidth  b,  as  shown  in  Fig.  3a.  The  power  spectral 
density  of  x^2(t)  will  then  be  as  shown  in  Fig.  3b*  [26],  The  multi- 
pliers remove  the  frequencies  centered  about  2fc.  Hence  the  signal 

has  spectral  density  shown  in  Fig.  3c. 

Thus,  each  term  x^(t)  Xj(t)  Eq.  (8)  contains  a d.c.  com- 
ponent (the  inpulse  at  w * 0 in  Fig.  3c)  and  a time-varying  component 
(the  power  contained  In  the  continuous  spectrum  from  0 to  b in 
Fig.  3c).  From  the  asynptotic  solution  to  Eq.  (8)  given  in  Eq.  (37), 
hoover,  it  is  clear  that  only  those  frequency  conponents  of 
x-j(t)  xj(t)  that  lie  below  the  cutoff  frequency  w0  will  have  an  effect 
on  the  solution  to  Eq.  (8).  We  may  safely  ignore  all  frequency  com- 
ponents of  x-j(t)  xj(t)  above  w0.  Returning  to  the  Ohio  State 
adaptive  array  as  an  example  (with  10  M-Hz  bandwidth),  we  see  that 
the  time-varying  portion  of  x-j(t)  xj(t)  will  have  baseband  frequency 
components  from  0 to  10  M-Hz.  The  feedback  loop  bandwidth  (w0)  is, 
however,  only  1 k-Hz.  Thus,  the  total  power  In  the  time-varying  part 
of  xn-  (t)  xj(t)  is  only  about  10  * of  the  d.c.  power.  The  part  from  0 to 
1 k-Hz  need  not  be  included  in  the  differential  equations  (8).  Re- 
ferring to  the  one-dimensional  differential  equation  of  Eq.  (20),  we 
see  that  the  situation  is  equivalent  to  having  8«  A.  In  this  case, 
we  may  safely  neglect  the  B term  in  constructing  the  solution  to 
Eq.  (20). 

*In  general  xj(t)  = Xj(t-tj),  where  tj  is  some  time  delay.  Thus, 
xi(t)  xi(t)  = Xi(t)  xj(t-Ti).  For  this  discussion,  we  assume  t,-  = 0, 
so  x.j  (t;  xj(t)  = x1,2(t).  The  value  of  tj  Is  not  inportant  for  the 
arguments  advanced  here,  which  are  only  qualitative. 
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Thus  we  will  solve  Eq.  (8)  by  including  only  the  constant  (d.c.) 
terms  of  the  products  x^(t)  Xj(t).  To  reiterate  the  arguments  above, 
the  validity  of  this  procedure  rests  on  two  observations.  First,  the 
solutions  to  Eq.  (8)  exhibit  a "cutoff  frequency"  effect.  Time- 
varying  terms  in  the  products  x^(t)  Xj ( t ) whose  frequencies  are  higher 
than  this  cutoff  do  not  affect  the  solution.  Second,  in  practical 
array  designs  for  communication  systems,  the  value  of  this  cutoff 
frequency  is  very  small  compared  to  the  RF  bandwidth.  As  a result, 
the  power  contained  in  the  time-varying  part  of  x^(t)  xj(t)  below  the 
cutoff  is  extremely  small  compared  to  the  power  in  the  constant  term. 
For  this  reason,  we  may  neglect  the  time-varying  portion  of  Xj(t)* 
xj(t)  entirely  in  solving  Eq.  (8). 

When  the  element  signals  are  random,  the  multiplier  output  is 
the  d.c.  portion  of  the  product  of  two  random  processes.  We  note  that 
this  quantity  is  the  same  as  the  infinite  time  average  of  the  product 
of  the  two  processes. 

Having  shown  that  only  the  constant  part  of  x^(t)  xj(t)  needs 
to  be  included  in  Eq.  (8),  we  further  note  that  the  constant  part 
may  be  obtained  from  the  infinite  time  average  of  x^(t)  xj(t).  More- 
over, when  all  x^t)  are  assumed  to  beergodic  processes,  the  time 
average  may  be  replaced  by  an  ensemble  average.  We  now  return  to 
the  discussion  of  Eq.  (8).  If  we  define  the  matrices 
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then  Eq.  (8)  can  be  written  in  matrix  form  as 


(42)  ^ + 2k*  W = 2kS. 

Let  us  consider  the  response  of  Eq.  (42).  When  there  is  more 
than  one  antenna  element,  the  system  of  equations  will,  in  general, 
be  coupled  in  ♦.  In  order  to  solve  this  system,  we  first  make  a 
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rotation  of  coordinates  Into  the  principal  axis  of  ♦.  Let 


] 

(43)  W = Rn 

1 

■j 

. where  R is  a 2N  x 2N  orthogonal  coordinate  rotation  matrix. 


represent  a new  system  of  coordinates  for  the  weights.  By  substi- 
tuting Eq.  (43)  into  Eq.  (42)  and  multiplying  on  the  left  by  R"1, 
Eq.  (42)  becomes 

(46)  ^-+  2k[R_1*R]n  = 2kR-1S  . 

If  R is  chosen  so  that  R-1$r  Is  diagonal. 


21 


(47)  R“'*R  = A = 


0 a2  o 


then  the  components  of  n lie  along  the  principle  axes  of  * and  the 
system  of  equations  is  uncoupled.  We  define 


(48)  P = R~  S = P 


and  Eq.  (41)  becomes  simply 


(49)  |n.  + 2kAn  = 2kP 


We  refer  to  the  components  of  n as  the  "normal  weights"  of  the  array. 

The  form  of  the  general  solution  to  Eq.  (49)  depends  on  the 
matrix  <J>.  Since  $is  real  and  symmetric,  its  eigenvalues  are  neces- 
sarily real.  Furthermore,  4>is  non-negative  definite  [2i] . Since 
none  of  the  eigenvalues  of  4>  can  be  negative,  the  solutions  to 
Eq.  (49)  will  not  contain  any  exponentially  growing  terms.  Further- 
more, none  of  the  eigenvalues  can  be  zero  when  there  is  element 
noise  in  the  arrqy  [21]  . By  "element  noise"  we  mean  random  noise  due 

J 

to  RF  components  behind  each  element  of  the  array.  This  type  of 
noise  is  incoherent  from  one  channel  to  the  next.  Element  noise 
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should  not  be  confused  with  a directional  noise  signal  received  by 
the  array;  such  noise  would  be  highly  correlated  between  elements. 

When  element  noise  is  present,  the  element  signals  (following 
the  quadrature  hybrids)  are  of  the  form 


(50)  x^t)  = ni(t)  + s . (t) 

where  n-j(t)  is  the  noise  component  and  s-j(t)  is  the  received  signal 
component  of  Xj(t)  (see  Figs.  2).  When  this  x^t)  is  substituted 
into  Eq.  (39),  «>is  found  to  be 


(51)  » = Vn 


a2 1 + <t> 

n s 


where  we  use  ^+n  to  denote  $ when  both  signal  and  noise  are  present 
and  when  only  signal  is  present.  on2  denotes  the  mean-square  value 
of  nj (t) : 

(52)  n-(t)  = 


(we  assume  all  n^t)  have  the  same  mean-square  value),  and  I denotes 
the  identity  matrix.  To  derive  Eq.  (52),  we  have  made  use  of  the 
assumption  that 


(53)  n.j(t)nj(t)  = 0 for  i^j. 
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- 


and 


(54)  ni(t)Sj(t)  = 0 for  all  i,j. 

Since  the  matrix  0n2I  1S  unaffected  by  a transformation  of 
the  type  R” 1 (o n2 1 ) R * °n2I,  the  same  orthonormal  matrix  which  diagon- 
alizes ^ will  also  diagonalize  $s+n.  Hence,  each  eigenvalue  of 
must  be  equal  to  on2  plus  the  corresponding  eigenvalue  of 
The  form  of  4>,  the  eigenvalues  and  the  rotation  matrix  will 
also  be  dependent  upon  the  number  of  signals  incident  on  the  array. 

If  we  assume  the  incoming  signals  to  be  uncorrelated  with  one  another 
and  the  element  signals  to  be  of  the  form 

n 

(55)  x^t)  = ni  (t)  + ^ sij(t) 

then  <^+n  will  be 

(56)  V„  ’ “S'  * *sj 

where  is  the  ^ when  only  signal  Sj(t)  is  incident  on  the  array. 
Since  the  rotation  matrix  R will  be  made  up  of  the  eigenvectors  of 
4>s+n,  its  determination  will  become  increasingly  difficult  as  the 

J 

number  of  signals  increases.  A method  of  determining  R and  A for  an 
arbitrary  nunber  of  signals  is  presented  in  Appendix  I.  We  will, 
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however,  in  the  next  section  present  a transformation  which  will 
allow  solutions  for  the  weights  to  be  determined  and  the  order  of  the 
system  reduced  by  a factor  of  2 with  an  arbitrary  nurrter  of  signals 
present. 


III.  THE  GENERAL  SOLUTION 

The  two-element  adaptive  array  behaves  in  a manner  similar  to 
higher  order  systems  in  that  it  has  pattern  flexibility.  It  is  there- 
fore capable  of  suppressing  directional  interference  signals  and  its 
anti -jam  performance  may  be  studied.  In  this  section  we  develop  the 
general  weight  solutions  for  a two-element  adaptive  array  assuming 
CW  signals  and  element  noise  are  present  at  the  array  input. 

The  system  of  equations  describing  the  response  of  a two- 
element  adaptive  array  is  from  Eq.  (42), 


Wj(t) 

w2(t) 

w3(t) 

w4(t) 


+ 2k  « 


Wj(t) 

w2(t) 

W3(t) 

W4(t) 


( s.<0  ^ 

s2(t)  I 

s3(t)  j 
Mt) 


When 

the  signals 

i n ci 

dent  on 

the 

shown 

that 

: ■twill 

be 

of  the 

form 

( $11 

$12 

$13 

$14 

(58) 

4>  = 

"$12 

$11 

” $ 1 4 

$1  3 

$31 

$32 

$33 

$34 

1 “ $ 32 

$31 

"$34 

$33 

It  has  been  demonstrated  by  Compton  [21]  that  when  $ exhibits  the  type 
of  symmetry  of  Eq . (58),  the  system  of  Eq.  (57)  can  be  reduced  to 


the  form 


1 , . d 1 

f wrjw2  ^ 

1 

j <59>  at 

^ W 3” JW4  j 

+2k 

♦ 1 1 + J ♦ 1 2 $1  3+J4>l4  W Wj-JW2 

♦ 3 1 + J <*> -»2  *3  3+jl’3it  J l Wa-jw,, 


= 2k 


(we  will  use  w-(t)  and  Wj  interchangably  throughout).  We  have  then 
formed  a set  of  differential  equations  in  terms  of  complex  weighting 
coefficients.  Recall  thatWj  and  are  the  in-phase  weights  and  w2 
and  w4  the  quadrature  weights  of  the  array.  The  quadrature  channel 
represents  a 90°  phase  delay  over  all  frequencies.  This  quarter 
cycle  delay  is  symbolized  mathematically  by  the  operator  -j. 

The  general  configuration  of  a two-element  adaptive  array  is 
shown  in  Fig.  4.  Let  us  suppose  these  are  n CW  signals  incident 
upon  this  array  from  angles  e,-  off  braodside.  We  will  assume  m of 
these  signals  to  be  desired  and  (n-m)  to  be  directional  CW  jammers. 
This  situation  is  illustrated  for  n=2,  nrl  in  Fig.  5. 

We  define 


(60)  Wj  = w i ~ jw2 


and 


(61)  w2  = w3-jw4 


27 


FEEDBACK 


ERROR  REFERENCE 


Fig.  4 — Two  element  array. 


DESIRED 

SIGNAL 


Fig.  5— -Signal  angles  in  space. 
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where  indicates  a complex  variable.  We  also  denote  the  colurm 
vector 


(62)  w 


\*2  J 


The  complex  correlation  matrix  Q,  will  be  defined  as 


(63)  $ = 

Vc  i>  J 

where 

(64)  A = , 

»n  + 

(65)  B = $i3  + j<f>i4 

(66)  C » 031  + j 4>32 

(67)  D = 4>3  3 + j 4>34  . 

Finally,  if  we  define  the  vector  S to  be 


(68)  S = 


( sr  js2  ^ 

Si  1 

Is  \) 

^S3-  js4  ) 
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then  Eq.  (57)  can  be  written  in  matrix  form  as 


, . H ^ 'V.  A.  a. 

(69)  w + 2kQ  w s 2kS  . 

Assuming  CW  signals  are  present,  the  received  signals  will  be 
of  the  form 

(70)  s.(t)  = • os.  • COS(w.jt) 

where 

o$^  = the  signal  amplitude 
and 

a>i  = the  signal  radian  frequency. 

Specifically,  at  the  second  antenna  (see  Fig.  5),  the  signal  is  given 
by 

n 

(71)  y2(t)  = l fz  • a • cos  ( <»> . t ) 

i=l  5i  1 

and  at  the  first  antenna, 

(72) y  (t)  l JT  • o . cos (w. t-o. ) 

1 i=l  si  1 i > 
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where  Is  the  phase  shift  between  antennas  due  to  the  propagation 


delay, 

(73)  a.-  = sin  e. 

1 Aa^j  1 


(L  is  the  element  spacing  and  xu1-  is  the  free-space  wavelength  at 
frequency  .) 

With  element  noise  present,  the  in-phase  and  quadrature 
element  signals  are 


(74)  xi(t)  = $i(t)  + m(t)  = l os. 

i=l  1 

(75)  x2 (t)  = s2(t)  + n2(t)  = ? o«. 

i=l  i 

(76)  x3(t)  = s 3(t)  + n3(t)  = J os. 

i=l  1 

(77)  x4(t)  = s4(t)  + n4(t)  = l as . 

i=l  1 


COS  (w^.  t-a^.  ) + n j (t) 
sin(u^t-o^ ) + n2(t) 
cos(u^t)  + n3(t) 
sinC^.t)  + n4(t). 


The  factor  was  included  in  Eqs.  (71)  and  (72)  to  make  the  in-phase 
and  quadrature  signals  in  Eqs.  (74)  - (77)  have  unit  amplitude. 
Substituting  Eqs.  (74)  - (77)  into  Eq.  (39)  yields  for$: 
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L 


* l 


\ l Oj  -coso.  \ l °l  sina\ 

2 1*1  si  1 2 i«l  st  \ , 

\ [ 

i J,  ”1,  ? j,  i ' 


? j,  % “s“<  0 ? j,’*,*0" 


We  observe  « exhibits  the  symmetry  of  Eq.  (58). 

Suppose  furthermore  that  the  reference  signal  is  given  by 

m 

(79)  R(t)  = A • l cos(u..t) 
i=l  1 

where  A is  some  non-zero  constant.  In  a practical  system  where  the 
reference  signal  is  derived  from  the  desired  part  of  the  array  output, 
A will  be  chosen  to  have  a value  compatible  with  the  equipment.  This 
fixed  operating  level  is  generally  achieved  by  inserting  a limiter 
somewhere  in  the  reference  signal  generating  network  | 15,25]  . In 
Eq.  (79),  the  reference  signal  is  coherent  with  the  first  m signals. 
If  Eqs . (74)  - (77)  and  Eq.  (79)  are  substituted  into  Eq.  (40),  we 
find 


r 

L 


L 

r 

• 

L 

l 

r 

w 


r 

L 


32 


(80)  S 


f*l\  , 

( ™ n \ 

s = 

1 1 

s2 

= A 

z cos(w^t)  2 cos(w.t-o.)  + m(t) 

i=l  i=l  si  1 1 

S3 

| 

m n 

i £ c°s(«it)  I o s1n(«it-oi)  + m(t) 

V i=l  i*l  i , 

1 

4 


Since  the  reference  signal  is  uncorrelated  with  the  directional 
jammers  and  element  noise,  we  obtain 

f .1  0 c cos  o,A 


(81 ) S=  | 


0e  cos  a. 
i=l  si  1 


" I °c  sin  a. 
i=l  si  1 


? 


v ’v  7 

Substituting  Eq.  (81)  into  Eq.  (68)  yields  for  5: 

M 


(82)  S = 


f Si-js2^ 


\ S 3“ JS 4 J 


( C 

Z °s  *e 
i=l  si 

m 

\k\  J 
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We  can  also  obtain  the  reduced  form  of  ♦.  From  Eqs.  (64) 


(67)  and  Eq.  (78)  we  obtain 


The  system  of  equations  describing  the  response  of  the  array 
may  now  be  determined.  Making  the  appropriate  substitutions  into 
Eq.  (69)  yields 


(84) 


[o2  +o2 

1=1  Si  ° 


n „ -J  a . 

Z °c  e 

i=l  si 


/ ? joA 

1 l a e ' 
i=l  si 


m 

V 1 °S 

V*’  1 j 


We  notice  that  the  equations  are  coupled  in  Q.  Solving  the  system 
will  require  a rotation  into  the  principal  axes  of  Q.  Had  the  order 
of  the  system  not  been  reduced,  we  would  have  a system  of  four 
coupled  differential  equations  to  solve.  The  coordinate  rotation  to 
uncouple  the  system  would  be  extremely  tedious.  Diagonalizing  the 
system  of  Eq.  (84)  will  be  obviously  much  simpler. 

We  begin  by  determining  the  eigenvalues  of  Q.  Since  0 is 
Hermitian,  its  characteristic  values  are  necessarily  real.  Let 


34 


J l ^ - * M a?  e _1 

2 i=l  si  n 2 i=!  si 

\ I °l  e''01  J-  [ o2  + a2  - X 

* i=l  si  2 i=l  si  n 


(85)  IQ-AI  | = 


- 0 . 


yields 


(I 

Hereafter  ^ will  be  represented  by  Expanding  Eq.  (85) 


(86)  [l  j°s.  +°n‘AJ‘i(|% 


1 lx  2 "j°‘.  r 2 n 

I U V e l °c  e = 0. 


v S . 

1 1 


Solving  the  above  quadratic  equation  for  x gives  for  the  eigenvalues 
of  Q: 

(87>  + * T 

i7j 

and 

(88)  X2  = o2+  1 Jo2  - I l o'(  +2  £ l o2  Oj  COS(a^  -a  • ) ] 
i i c \i  si  1 J i j 1 J / 


From  Eqs.  (87)  and  (88)  we  can  now  determine  the  eigenvectors 
of  Q.  Let  Ej  denote  the  eigenvector  associated  with  x^ . We  require 

i 

(89)  (Q  - *.!)(£.)  = 0. 


] 


Substituting  Eqs.  (83)  and  (87)  into  Eq.  (89)  yields 


(90) 


V 621 

We  then  obtain 
(91)  }Tll  + 


and 


(92)  G*en 


where 


(93)  N - \ 


(94)  = \ 


and  superscript 
notice 


ll  o]  o?  COS(a..o.) 

i j i j J J 
i^j 


\\/2 


f ( 


\ 

7 


+ 2 l JV  02  cos( 
i j i i 
i^j 


two  equations  of  the  form 


V\/ 

G e2 i = 0 


= 0 


°o  al  COS (a. -a.) 

i j 1 J 


"*"  denotes  the  complex  conjugate.  However,  if  we 


p ! 


] 

] 

] 

] 


] 


1 

J 


] 


] 


] 


•V,  'lA. 


(95)  F2  = BB* 


then  the  eigenvectors  can  be  chosen  to  be 


(96)  E 


\l/2 


A 


i?j 

-oa,- 


J 


We  hereafter  denote  £ £ by  £ J. 

i=l  j=l  ij'j 

Similar  calculations  show  the  second  eigenvector  to  be 

2 y %% cos(ai -“j1)172  i 


(97)  E,  = 


1 "Ja-i 

1 V 2 „ 1 

* * K * 


If  from  the  eigenvectors  we  form  the  matrix 
(98)  Z = (Ej,E2) 
then 


1%  % 

(99)  Z Q Z = A 
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where  A is  a diagonal  matrix  composed  of  the  eigenvalues  of  Q and 
superscript  1 denotes  the  inverse  matrix.  If  from  2 we  construct 
a unitary  matrix  R given  by 

(100)  R = 1 

ik 

then 

, i 

(101)  R = R 

where  "+"  represents  transpose  conjugation.  We  require 

(102)  R'R  = I 


where  I is  again  the  identity  matrix.  From  Eqs.  (96)  - (98)  and 
Eq.  (100),  we  obtain  for  R: 


f r* 

Performing  the  operation  indicated  in  Eq.  (102)  yields  for  the  1 : 

normalization  constant  k:  r 

r 

i 


[ 
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and  using  the  correspondence  of  Eq.  (99),  we  can  rewrite  Eq.  (106) 
conveniently  as 


(109)  {!?■  r + 2kA?  =2kV. 


dt 


The  elements  of  fwill  be  referred  to  as  the  "unitary  weights"  of  the 
array.  This  is  analogous  to  the  "normal  weights"  of  the  array 
represented  by  n in  Eq.  (49),  when  all  the  elements  are  real. 

Recall  that  A is  a real  diagonal  matrix  and  is  composed  of  the 
eigenvalues  of  Q.  Specifically  from  Eqs.  (87),  (88)  and  (99),  we 
find 
(110) 


A = 


( ? K*0"*  * 2 y^ycos(vv) 


l 

\2 


V 


5"  H +2  l [°c  °c  cos  (a.  -a . ) )2  j 

2 i si  n 2 Vi  si  y si  sj  1 yJ  J 


Since  r and  V represent  coluim  vectors,  the  equations  are  no  longer 
coupled.  If  we  now  carry  out  the  indicated  operation  of  Eq.  (108),  we 
obtain  for  V 


(111) 


fl  o e^V  7-i-=1  $i  ^ $i  ~ 

r1  1 IK  + *U°L*l. 


i.  ° * i 


Ja. 


V = R S = 


2^2] 


cos (a. -a.) 

i«  Ji  °j  1 J 


0/2 


m fa  \ a*  e^ 

l a ja,-  1*1  si  1-1  Si 


V=1  si e 


(f  °*i  + 2 °S.°S.  cos(ar°j> 


11/2 
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We  are  now  in  a position  to  obtain  the  equations  describing  the 
unitary  weights  of  the  array.  Making  the  appropriate  substitutions 
into  Eq.  (109)  yields  a set  of  separate  equations  which  may  be 
written  as 
(112)  . 

1 

Ht  Yi<t)+2k[an+  \ ?as.+  +2^  £as  as  cos(Oi-o.)  fjy^t)  = 

i 5i  ‘Vi  si  i/j  si  sj  1 J ' 


m n ja. 
m ja.  k .Z  °s .1  °s.e  1 

k .A.  I o.  e 1 + — '-=1  11  = 1 i 


i=l  si 


and 

(113) 


(la s.  + 2 l la s °c  COs(a  -o j)2 

M si  i/j  si  sj  1 J / 


3&<‘>  + 2kK*  ? H.  - H K * 2 l H °|  cos(a.-a  )¥]?2(,)  . 

'it  i£j  i j J J 


ffl  J°4 
k -A-  l a e 
i=l  si 


m n ja. 

k l o l o2  e 1 
i=l  Si  i=l  si 


(I  °s  + l I°f  °s  COS  (a  -a  .)¥ 

m si  i7j  si  sj  1 J r 
The  solution  of  Eq.  (57)  has  essentially  been  reduced  to  solving  two 
first-order  differential  equations  of  the  form 


(114)  3t  Yi(t)  + Vi(t)  = kV 
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; 


Equations  (112)  and  (113)  are  then  easily  solved,  yielding  as 
solutions: 


(115) 


r,(t)  = kje 


/ m n ja. 

A I ? 6 

+ 2 as  e + 

2-  2 i=l  si  i 

J v l r i i 


(l°S  +2^  1°S  as  COS(oi-ai)V 

Vi  5i  i^j  5i  5j  3 > 


(116) 


•V  -2k  A V A m 

y2(t)  =k2e  +~^~  l a 

2 A?.2-  2 i=l  si 


m ja. 

2 °c  *e 


m n ja. 

I °s  l al  e 

i=l  5i  i=l  5i 


+?l  ^2  COS (a. -a.] 

Vi  si  Si  Sj  1 J 


where  kj  and  k2  are  the  constants  dependent  upon  the  initial  weight 
values  and  Atand  a2  are  the  eigenvalues  of  Q given  in  Eqs.  (87)  and 
(88).  It  is  apparent  that  the  complex  weights  and  W2  may  be  ob“ 
tained  from  the  unitary  weights  Yj  and  y2  by  the  simple  transformation 

'll  *\» 

(117)  U = R r . 


Then  from  Eq.  (105),  we  find 


•V  \ + \ 


(118)  w = — y 


The  following  relations  are  easily  verified: 


(122)  *2-A!  = “ U °s  + 2 H °s  cos(a.-a.) 


i i 


(123)  A1+X2  = 2 02  + I o2 


l 1 


and 


(124)  4-Xj-X  = 4o2(02  + £0*  ) + 2 £ l o|  o|  (1  - COS(a.-a,)) 

i i i j i j 


Using  the  foregoing  results,  Eq.  (120)  readily  reduces  to 


-v,  ■■  -2kX  t -2kX  t 

(125)  w (t)  = — [ k e + k e ] 

1 2 


+ A 


?!  a 


JiV  'jUr^^J'lJisyrsV  i 


4an(on  + laV  + 2 ' l H as.(1“cos(ai'aj)) 

1 J * J 


1 1 


i 


j 

1 

'J 


J 


+ A 


^an(°n  + 

n n 


I°l  ) + 2 I Ia|  a|  (1 -cos (a. -a.) 
i i i j i j J 


From  Eqs.  (60)  and  (61)  it  is  seen  that  the  real  array  weights 
can  simply  be  found  from  the  complex  weighting  coefficients  by  means 
of  the  following  relations: 

(127)  Wl(t)  = Real  (wx(t)) 

(128)  w2(t)  = -Imag  (Wj(t)) 

(129)  w3(t)  = Real  (w2(t)) 

(130)  w4(t)  =-Imag(w2(t)). 

To  summarize,  then,  we  have  developed  expressions  for  the 
weighting  coefficients  of  a two-element  adaptive  array  with  an 
arbitrary  number  of  jamming  and  desired  CW  signals  incident  upon  it. 
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Let  us  examine  the  temporal  response  of  the  weights.  From 
Eq.  (125)  and  (126)  the  transient  response  of  the  complex  weights  are 


^ , -2kx.t  -2kX  t 

(131)  »ltr.U)  - — ( kj  • e + k2  • e ) 


and 


(132)  wnr_(t) 


2l  H COS(<y 

i^j  i J 


aj) 


-2kx  t -2kX  t 
kje  -k2e 


We  notice  the  transient  response  contains  decaying  exponentials.  The 

rate  of  decay  of  these  functions  is  proportional  to  the  loop  gain 

% 

constant  k and  the  eigenvalues  of  Q,  \l  and  \2.  Recall  that  these 
eigenvalues  are  functions  of  the  signal  powers  oSi2,  the  noise  power 
on2  and  the  electrical  phase  shift  between  antenna  elements  . We 
then  observe  that  as  the  number  of  signals  increases,  the  eigenvalues 
also  increase  (See  Eqs.  (92)  and  (93).)  The  array  response  will  then 
become  faster  as  the  number  of  signals  increase.  Also,  a larger 
element  noise  power  a n2  and  loop  gain  constant  k results  in  a faster 
array  response.  The  array  weights  then  converge  to  their  final 
solution  in  a shorter  period  of  time  with  many  signals  present  than 
with  few.  Unfortunately,  as  it  will  be  demonstrated,  the  output  SNR 
degrades  as  the  number  of  signals  increases. 
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IV.  THE  PERFECTLY  CONSTRAINED  ARRAY 


The  results  of  the  previous  section  may  be  used  to  study  the 
performance  of  a two-element  array  when  there  are  two  signals  present 
When  an  N-element  array  has  N signals  incident  upon  it,  the  array  is 
said  to  be  "perfectly  constrained."  It  will  be  shown  that  with 
element  noise  present,  the  SNR  at  the  array  output  is  degraded.  Not 
only  is  the  desired  signal  amplitude  below  that  of  the  reference 
signal  1 16]  but  the  interference  signal  is  also  present  at  the  array 
output. 


Let  us  suppose  two  CW  signals  are  incident  upon  the  array 
(see  Fig.  5).  The  signal  arriving  at  angle  0j  off  broadside  is 
chosen  to  be  the  desired  signal  and  the  other  is  a directional  inter- 
ference signal.  From  Eq.  (125)  - (130)  the  steady-state  weights  are 
found  to  be 


(133)  Wl(») 


2 o2(a2+o2  +o2  ) + 2o 
n'  n si  s2' 


2 • cos  a?) 

5 2 

5lOS2^'C0S(°1-a2)) 


(134)  „ (.)  = ~ ■ si"  °l> 

4an(an+as1+ai2)  + 


A-oSi(2o^+o2?»  (l-COS(ai-a2))) 

4an  (%+as  1+0s2)  +2as  !as2  (1  -cos  ( °i-°2 } } 


and 
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-l 


I 


(136)  w (-)  = 


a’°s1°s2  * s*n(°i-«2) 


4on(on+os1+os2)+2os1°s2(1-c0s(ai'a2)) 


The  output  of  a two-element  array  with  n CW  signals  Incident 
upon  it  is 

n 

(137)  s(t)  = o$  ((w^  + wj  S1n  ai  + w2  cos  °i) 


2 p 

+ (w  + w cos  a.  - w sin  a. ) ) • cos(w. t 

w i 1 Z I 


- tan 


-if  W4+Wl  s^n  ai  + wi  cos  ai 

^ w3+Wj  cos  a7  - w2  sin  jj 


Substituting  Eqs.  (133)  - (136)  into  Eq.  (137),  we  find  that 
the  desired  part  of  the  array  output  is 


A'O^  • (4a^+2a^(l-COS(a1-a2)))*cos(Wlt) 


(138)  Sd(t)  = 


4on(an+°s  1+as2)  + 2a^o^(l-cos(ai-a2)) 


2 2 


We  notice  that  because  element  noise  is  present,  the  desired  signal 
output  does  not  match  that  of  the  reference  signal  given  in  Eq.  (50). 
Similar  results  were  obtained  by  Compton  [16]  . However,  this  is  not 
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the  only  effect  of  the  element  noise.  If  we  calculate  the  undesired 
part  of  the  array  output,  we  find 


J 


] 

] 

] 

I 

J 

] 

] 

1 

] 

] 

1 

] 

] 

] 

] 

] 

] 


(139) 
S,(t)  = 


2 

A’o.  o',  o • cos 
Si  $2  n 


a2~a2 


4on  K + °l1  + as2}  + 2aS!  + °s2(1  " C0S  (“ra2^ 


f 

- ll 

f sin(a2-ai) 

> 

COS 

w2t-tan  I 

\ 

\ 

l l+COS(a2-cti)  j 

/ 

We  see  that  the  jamming  signal  is  also  present  at  the  array 
output.  By  changing  the  value  of  the  weights,  the  element  noise 
causes  the  spatial  null  formed  by  the  array  weights  to  no  longer 
exactly  be  in  the  jamming  signal  direction.  For  example,  in  Figs.  6 
and  7 we  see  the  array  pattern  with  one  jamming  and  one  desired 
signal  present.  In  both  cases  the  jamming  and  desired  signals  have 
equal  power.  In  Fig.  6,  the  jamming  signal-to-element  noise  power 
ratio  (ctj 2/crn2)  is  20  dB.  The  element  noise  has  little  effect  on  the 
array  pattern.  The  jamming  signal  is  well  within  the  null.  However, 
as  the  element  noise  power  is  increased,  the  null  moves  farther  away 
from  the  jamming  signal  direction.  For  example,  in  Fig.  7 we  see  the 
pattern  when  aj 2/on2  = 6 dB  (the  desired  signal  power  and  interference 
power  are  unchanged).  We  notice  the  jamming  signal  power  is  now  only 
18  dB  below  the  desired  signal  power.  This  may  be  of  serious  con- 
sequence when  both  the  jamming  and  desired  signals  are  low-power 
signals  (power  levels  near  the  element  noise  power).  We  notice  from 
Eq.  (139) 
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(140)  Tim  S-(t)  = 0 


In  other  words,  as  the  element  noise  power  becomes  small,  the  jamming 
signal  is  no  longer  present  at  the  array  output. 

The  situation  may  also  arise  that  both  incoming  signals  are 
desired  |15|.  Similar  calculations  show  that  the  desired  signal 
amplitudes  are  again  below  the  reference  signal  amplitude  when  element 
noise  is  present.  The  phase  of  the  desired  signals  is  correct. 

We  have  seen,  then,  that  the  array  makes  a compromise  between 
the  contribution  to  e2(t)  due  to  the  element  noise  and  those  due  to 
the  signals.  Since  the  element  noise  is  uncorrelated  between 
channels,  the  array  is  unable  to  remove  its  contribution  at  the 
output. 
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V.  THE  OVER- CONSTRAINED  ARRAY 

The  results  of  section  III  may  also  be  used  to  study  the 

t 

performance  of  a two-element  array  when  more  than  two  signals  are 
present.  When  an  N-element  array  has  more  than  N signals  incident 
upon  it,  the  array  is  said  to  be  "over-constrained." 

Let  us  assume  a two-element  array  has  three  CW  signals 
incident  upon  it.  In  other  words,  the  antenna  signals  are  from 
Eqs.  (71)  and  (72): 

(141)  Yl(t)  = 2{aSi*cos(w1t-a1)+aS2*  COs(w2t-a2)  +0^*  COs(w3t-a3)} 
and 

(142)  y2(t)  = ^^°Sl  * cos (w jt ) + • cos(w2t)  + o$^  • cos (w3t) } . 

First,  consider  the  response  with  no  element  noise  present. 

We  assume  the  signal  at  frequency  Wl  to  be  desired  and  the  others  to 
be  CW  interference  signals.  From  Eqs.  (125)  - (130)  we  determine  the 
steady-state  weights  to  be: 

A-°srK2(c0S  ai  ■ cos  a2 ) + (cos  aj-  cos  a3) 

(143)  w ! ( °°)  = — _i_ 

2(aSi'OS2(l'C0S(“l-a2))+Os2OS3(1'C0S(a2‘a3))+os1as3(1"C0S(ai-O3)))  i 
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! 


A*a  {a2  (sin  az  - sin  oti)  + a2  (sin  03  - sin  a;)) 

Si  S2  b 3 

2(o*  o*  (l-cOS(a ,~a))+a2  a2  (l-COS(a  -a  ))+a^  0 2 (1  -COS (a  -a  .) ) ) 
SlS2  IS  5 2 ^ 3 51  s3 


A-a  • (a2  (1 -cos (a . -a  ) ) + a2  (1 -COS  (o^-cij) ) ) 

s l s2 _3 

2 (of  (l-COs(a,-a  ))+a*  o (l-COS(a.-ct3))+a*  a2  (l-cosfaj-a  ))) 
S 1 S 2 S2  3 1 3 


(146)  w4(“) 


A-o  (<r  • sin(a  -a  ) + a • sin(o  -a  )) 
S1  s2  b3 


2(o2  (l-cos(ara2))+Oj  a2  (l-cos(a2-a3))+a*  0^  (1  -C0S (a^a 3) ) ) ’ 
12  _ , 13 


We  can  substitute  Eqs.  (143)  - (146)  into  Eq.  (137)  and 
determine  the  array  output.  For  the  desired  signal,  we  have 


(147)  sd(t)  = 


A‘as  (°S  ^ ’ C0S(ara2))  + °S  0 “COS  (a^-ctj)  ))  'COS  (w  jt) 


1 2 


o 2 -o2  (1 -COS (a  -a „))+o2  a2  (1 -C0S(a  -a  ) )+a*  a2  (l-COS(o  -a  ))  . 

si  s2  2 2 3 1 3 


From  Eq.  (147)  we  notice  that  even  though  there  is  no  element 
noise,  the  amplitude  of  the  desired  part  of  the  array  output  does  not 
match  that  of  the  reference  signal  (recall  the  reference  signal  ampli- 
tude is  A).  We  do  notice  that  the  desired  signal  phase  is  correct, 
however. 
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A similar  calculation  for  the  interference  signals  at  the 
array  output  may  be  made.  When  this  is  done,  we  find  that  both 
interference  signals  are  present  at  the  array  output.  Performing 
these  calculations  and  using  Eq.  (147),  we  determine  the  signal-to- 
jam  (interference)  ratio  (SQ/J0)  at  the  array  output  to  be 
(148) 


4'as1^s?(1'cos^ara2^+as3(1"cos(ai"a3^)2 

4(a£  +o£  )+2a£  •cos(a1+a2-2a3)+2a£  cos(a1+a3-2a2)+4(aVo^  )cos(a2-a3) 


VJo  = 


*2  °3  3 3 

„/  4 . 4 , 


s2  S3 


2(os3+as2)cos<ara3) 


With  two  signals  incident  on  a two-element  array,  the  signal 
to-jam  ratio  (with  an2  = 0)  may  be  shown  from  Eq.  (138)  and  (139)  to 
be  infinite.  The  additional  jamming  signal  has  caused  a degradation 
in  S,/J0. 

Suppose  we  now  consider  the  first  two  signals  at  the  input  to 
be  desired  signals.  By  again  calculating  the  array  output,  we  find 
that  not  only  are  the  desired  signal  amplitudes  degraded  and  the 
jamming  signal  present  at  the  array  output,  but  the  desired  signal 
phase  no  longer  matches  the  reference  signal  phase.  In  an  adaptive 
array  for  use  with  more  than  one  desired  signal,  phase  control  of  the 
desired  signals  may  not  be  possible.  For  example,  the  four-element 
array  described  in  [ 15 ] for  use  with  sensors  as  many  as  four  desired 
signals  may  be  present  at  the  array  input. 
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It  has  been  demonstrated  in  this  section  that  when  an  adaptive 
array  is  over-constrained,  the  array  performance  is  degraded.  We 
have  found  that  the  desired  signal  amplitude  (and  possibly  the  phase) 
no  longer  matches  the  reference  signal.  The  interference  is  also 
present  at  the  output.  With  element  noise,  the  performance  will  be 
further  degraded  although  this  degradation  will  be  small  when  the 
element  noise  power  is  much  smaller  than  +-he  signal  powers. 


VI.  ARRAY  PERFORMANCE  WITH  A CONTINUOUS  INTERFERENCE  SECTOR 

The  results  of  section  VI  for  the  array  performance  with 
discrete  signals  will  now  be  extended  to  the  case  of  a continuous 
sector  of  radiation.  It  has  been  shown  in  section  II  that  when  the 
signals  incident  upon  the  array  are  such  that  their  bandwidths  are 
greater  than  the  product  of  the  loop  gain  constant  (k)  and  the  signal 
powers,  the  time-varying  components  may  be  neglected  in  the  differ- 
ential equations  describing  the  weights.  In  this  section  we  assume  a 
continuous  sector  of  spatial  radiation  is  incident  upon  a two-element 
array.  Moreover,  we  assume  the  frequency  spectrum  of  this  radiation 
the  same  from  all  angles  and  the  bandwidth  of  this  radiation  is  wide 
enough  that  we  may  neglect  all  but  the  d.c.  components  in  the  weight 
equations.  We  also  assume  the  phase  shift  between  antennas  due  to 
propagation  delay  will  be  nearly  constant  over  the  entire  signal 
bandwidth.*  Then 

( 149)  Aw.  -v  for  all  i . 

With  this  approximation  we  may  proceed  directly  from  the  values  of 
Q calculated  in  Eq.  (83).  Recall,  the  terms  of  0 were  of  the  form 

(i5°,  ',=Hvs,n  e,k  • 

where  o$l-2  is  the  power  of  each  discrete  signal.  We  define  D(  Qj ) to 
*We  have  neglected  the  effect  of  the  wave  polarization. 


57 


I 
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be  a spatial  power  density  defined  over  e < e<  e2.  Also,  define 
a to  be  a partition  of  the  interval)  and  ||a||  to  be  the  "norm" 

of  this  partition.  Then  let  us  suppose  we  make  the  number  of  signals 
go  to  infinity  while  keeping  the  total  power  constant;  i.e.,  as 
1 1 A 1 1 * 0 as  n + 

lim  l f sin(i  Ae) ) • D(i  Ae)  • i Ae 
I 1 A | |->0  1 ^ 1 ' 

remains  finite.  The  sumnation  may  then  be  replaced  by  a Riemann 
integral.  In  other  words,  we  have  replaced  the  model  of  many  discrete 
signals  of  a relatively  small  difference  in  frequency  (Aw)  with 
respect  to  the  carrier  by  a continuous  spatial  distribution  of 
radiation. 

Practically  this  situation  arises  when  a waveform  irradiates  a 
rough  surface  (e.g.,  surface  of  the  earth,  ocean,  ionospheric  skip, 
etc.)  and  the  reflected  energy  appears  to  be  spread  in  angle  with  no 
specific  source  of  origin.  Another  example  would  also  be  sky  noise. 

With  this  in  mind  we  write  Eq.  (150)  as 

(151 ) q = j f|2-^sine  D(e  )d  e 

Specifically, from  Eq.  (151)  and  the  relations  of  Eq.  (64)  - (67), 
the  components  of  Q now  become: 

' -v 


. 
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(152)  B = D ( 0 ) e “ 


. 2ttL 

j j — sin  e 


'0  f 2 

(153)  C = | D(e)e 


-j  ^ sin  e 


(154)  A = D = j D(e)de  . 


The  expression  for  Q is  then 


(155)  Q = 


j sine 


D(e)  . 2ttL 

-j  -jp  sine 


Define  the  relations: 


(156)  u = sin  e , e = sin-1  (u)  , 


( 157)  dp  = cos  e • de,  de  = ^ ^ — 

cos  e i ■ 

Vl-P2 


1 


A 


058)  e 


- 2ttL 


059)  f(p)  = ^m_J 
Ji  - p2 


Using  the  previous  definitions,  Eqs.  (152)  - (154)  become 


(160)  B(-b)  = j f(p)ejSlJ  dp 


„ f 2 

(161)  C(b)  * J f(u)e'Jew  dp 


(162)  A = 0 = p f(u)dv  = B(0) 


where  B(e)  is  the  Fourier  transform  of  f(u). 
veniently  be  written  as 


The  matrix  Q may  con- 


(163)  Q = 


B(0)  B(-6) 


B(B)  B(0) 
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It  is  worth  noting  here  that  if  the  incident  signal  is  arriving 
from  a discrete  angle  e^,  then  we  may  represent  this  case  as 


(164)  D(e)  = o2  6(e-ei) . 


Substituting  Eq.  (164)  into  Eqs.  (160)  - (162)  results  in 


j ^=-  sin  G 


(165)  B = j o2%  e u 
°i  1 


j sin  0, 


= 2 Q w 


•6(e-e.)de  = e 
si 


j ^ sin  e 

y u\ 


( 1 66 ) C = o'  e w • <5(0-0. )de  = e 
^ * s • 


-j  sin  8, 


(167)  A = D = a?s  6(0-0.)de  = o* 


Then  by  allowing  discrete  signals  to  be  represented  by  Dirac  Delta 
functions,  the  continuous  channel  representation  readily  reduces  to 
the  discrete  case  previously  investigated. 

If  we  assume  the  impinging  radiation  sector  to  be  a jamming 
source,  then  a desired  signal  source  must  also  be  present  in  order 
for  the  weights  to  have  non-zero  solutions.  We  define  • 
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(168)  B(e)  = (3)  + Gd(B) 


where  B j ( b ) is  due  to  the  jammer  sector  of  radiation  and  Bd( a)  to 
the  desired  sector.  Since  the  interference  and  desired  signals  are 
uncorrelated,  we  may  add  the  contributions  of  the  jamming  and  desired 
signals  in  B(b).  Using  Eq.  ( 1 68)  we  are  now  able  to  write  Q as 


(169)  Qd+j  * Qd  ♦ Qj  - 


( 


B(0) 

b(-b) 

b(b) 

B (0 ) 

If  we  now  assume  the  desired  signal  to  be  a CW  signal  at 
frequency  u)d  located  at  a discrete  spatial  angle  ed,  then 


fa2 

1 s 


(170)  Q.  = 
d 


• e 
sd 


It  e'j6"d 


J 

where  from  Eq.  ( 1 56)  we  have 


( 1 71 ) ud  = sin  9d 

and  oSd2  is  the  desired  signal  power.  The  reference  signal  R(t)  will 
again  be  of  the  form  of  the  desired  signal,  then 


1 


i 
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(l72)  R(t)  = A • cos 


From  Eq.  (80)  the  vector  S is 


(173)  S = ~ 


o • e 


Figure  8 shows  a two-element  array  with  a sector  of  jamming 
radiation  and  a single  desired  signal. 


DESIRED 
SIGNAL  A0- 


JAMMING 
.RADIATION 
\ SECTOR 


Fig.  8— Array  signals. 
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We  notice  that  if  f(u)  in  Eq.  (156)  is  real,  then 

(174)  B(-e)  = B*(fl) 

and  if  f(u)  is  purely  imaginary 

(175)  B(-B)  = -B*(8)  • 


By  substituting  Eqs.  (170)  and  (173)  into  Eq.  (69),  we  find 
that  the  complex  weight  response  is  determined  by 


1 

^B(O) 

B ( -6 ) A 

(176) 

1 ^ 

+ 2k 

yB(8) 

B (0 ) ) 

= K-A-a 


f^\ 


V ) 


The  method  of  solution  is  the  same  as  that  of  the  previous  sections. 
We  first  determine  the  eigenvalues  of  Q to  be 


A A 


(177)  A = 


V 


Ax  0 
0 A 


1 

^ B(0)  + (B(8)B(-8))2 

1 

2 


V 


B(o)-(B(e)B(-e)r 


J 


Since  f(u)  is  a power  spatial  density,  then  it  is  necessary 
real  and  Q is  Hermitian  with  real  eigenvalues.  Moreover,  we  note 
from  Eq.  (17^  that 
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1 1 

(178)  (B(e)B(-B))2  = (B(e)B*(f3))2  = |B(B)|  . 


Using  the  results  of  Eqs.  (168),  (177)  and  (178),  the 
coordinate  rotation  matrix  becomes 


(179)  R = 


2 iBjJLL 

\ Ib(b)  I 


Z-BM. 
|B(e) | 


Again  we  transform  the  coordinates  into  the  principal  axes 
of  Q.  The  expressions  describing  the  unitary  weights  are  obtained  by 
eval  uating 

( 180)  \t  = R+S 


(181)  r = R+ft. 


Substituting  Eqs.  ( 177),  (180)  and  (181)  into  Eq.  (69)  yields 


d_  ^{t)  B(0)+|B(6)|  0 W'jU)  + ^d 

dt  y2(t)  +2k  0 B(0)-|B(3)| fc2(t)  J~z 


K-A-o.  • ft*"*  ♦ Sll) 


d . lie 
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the  system  of  Eq.  (182)  is  easily  solved  yielding  as  solutions 


(183)  y,  (t)  = k.e'2™0^  B^l  * + — & — 

2-/2(B(0)  + I B ( 3 ) | ) 


(184)  Y,(t)  = 


A*  a e 
sd  ' 


d B*(8 


. -2k[B(0)-|B(e)|]t+ ^ 

2e  2.j2[B(0)-|B(e)l] 


where  kj  and  k2  are  the  constants  dependent  upon  the  initial  condi- 
tions of  the  system.  From  the  inverse  transformation  given  by 
Eq.  (116),  the  complex  array  weights  may  be  obtained  from  the  above 
expressions.  Again 


O/  'V\, 

(185)  w = Rr 


Making  the  necessary  substitutions  we  find  the  expressions 
for  the  complex  weight  response  are 


1 f -2kXjt  -2kX2t 


+ k2e 


A'o  [B(0)e  a-B*(0)] 

sd 

4 [B2(0)  - J B ( 3 ) | 2 3 
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and 


(187)  w?(t)*^- 


2|B(e)| 


kje 


■2kAjt 


- k2e 


2kA2t 


+ 


A*o  (B(0)-B(8)eJ6yd) 
sd 

4(B2(0)  - 1 B ( 0 ) |2) 


These  specific  general  expressions  are  given  as  a function  of  B(b), 
the  Fourier  transform  of  f(u).  f(u)  is  directly  related  to  the 
spatial  power  density  (Eq.  (156)).  The  array  weights  given  in 
Eqs.  (186)  and  (187)  then  may  be  found  from  the  spatial  power  density 
D(  e)  and  the  resulting  adaptive  array  performance  may  then  be  examined 
analytically. 

In  the  next  section  we  treat  several  specific  examples  of 
continuous  jamming  radiation.  For  convenience  we  have  used  the  results 
for  many  discrete  signals  to  represent  this  continuous  sector. 


,1 


67 


rl 

l I 
. i 


VII.  RESULTS 


In  the  previous  sections  the  values  for  the  array  weights  of  a 
two-element  adaptive  array  were  calculated.  For  the  conditions  of 
multiple  signal  interference  or  a continuous  sector  of  interference, 
the  steady-state  array  weights  may  be  found.  Digital  computer 
nrograms  have  been  written  to  carry  out  these  computations.  The 
programs  calculate  not  only  the  final  weights,  but  also  the  resultant 
array  pattern  and  the  signal-to-noise  ratio  SNR  at  the  array  output. 
The  programs  are  given  in  Appendix  II.  In  this  section  we  discuss 
the  results  of  these  computations. 

For  a continuous  interference  sector  it  is  found  that  when 
the  desired  signal  is  separated  from  the  interference  by  more  than 
7 degrees,  the  SNR  degrades  as  the  interference  sector  becomes  larger. 
.'Jhen  the  separation  between  the  desired  signal  becomes  small  (<  7°), 
the  SNR  surprisingly  increases  as  the  interference  sector  becomes 
wider.  This  result  is  especially  interesting  and  will  be  discussed 
again  later. 

Figures9  through  15  show  some  typical  plots  of  the  final  array 
patterns.  In  these  plots  the  desired  signal  amplitude  (aSl)  is  1. 

The  total  interference  power  is  then  lOdB  greater  than  the  desired 
signal  power  at  each  element.  The  element  noise  amplitude  is  .01. 

From  Eqs.  (125)  and  (126)  for  the  final  weight  values,  the  antenna 
patterns  are  generated.  The  element  spacing  is  chosen  as  A/2  in  each 
case . 
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In  Fig.  9 the  desired  signal  arrives  at  -30°  off  broadside. 

The  interference  sector  is  from  0°  to  1°.  We  notice  most  of  the 
interference  is  well  within  the  null.  In  Fig.  10  the  interference  is 
increased  to  5°  and  in  Fig.  11  to  15°.  (The  interference  power 
remains  the  same  in  each  case.)  We  notice  the  center  of  the  null 
remains  at  the  center  of  the  interference  in  each  figure. 

In  Fig.  12  the  interference  sector  remains  0°  to  15°  while  the 
desired  signal  is  now  at  -15°.  We  notice  no  pattern  change.  In 
Fig.  13  the  desired  signal  arrives  at  -5°.  A close  observation  shows 
that  the  null  has  moved  slightly  to  the  right  and  the  pattern  remains 
beamed  upon  the  desired  signal.  However,  as  the  interference  sector 
is  decreased  to  5°,  the  desired  signal  suffers  an  attenuation  of  5 dB 
as  may  be  seen  by  comparing  Fig.  13  with  Fig.  14.  A further  attenua- 
tion of  2 dB  is  also  evident  in  Fig.  15  when  the  interference  sector 
is  decreased  to  1°. 

Figures  16  through  19  show  the  signal-to-noise  ratio  SNR  at  the 
array  output  as  a function  of  the  interference  sector  width  for  various 
desired  signal  positions.  In  these  figures,  the  interference  sector 
varies  in  width  from  0°  to  ANG°.  In  Fig.  16  the  desired  signal  varies 
from  -60°  to  -5°  and  in  Fig.  17,  -10°  to  3°.  When  the  desired  signal 
is  within  7°  of  the  interference  sector,  we  notice  an  increase  in  the 
SNR  as  the  interference  angle  increases.  The  SNR  increases  because, 
as  we  have  seen  in  Figs.  13  and  15,  the  array  causes  the  null  to  turn 
away  from  the  desired  signal  as  the  interference  sector  increases. 
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In  Fig.  18  we  notice  the  change  in  SNR  is  small  beyond 

(AN6=)  5°. 

In  Fig.  19  we  notice  that  beyond  15°  there  is  no  change  in  the 
SNR,  regardless  of  the  desired  signal  angle.  For  a jamming  sector 
beyond  15°,  an  SNR  degradation  of  8 dB  is  apparent  in  Fig.  19  when 
the  desired  signal  is  at  -60°.  We,  however,  notice  an  SNR  improve- 
nent  of  4 dB  when  the  radiation  sector  is  increased  from  .1°  to  30° 
for  a desired  signal  arriving  at  -5°. 
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VIII.  CONCLUSIONS 


Tiiis  report  has  presented  a study  of  the  behavior  of  a two-element 
adaptive  array  with  multiple  discrete  signals  and  a conti uous  inter- 
ference sector.  It  has  been  determined  mathematically  that  when  the 
signal  bandwidths  are  large  compared  to  the  feedback  loop  bandwidth, 
only  the  d.c.  terms  need  be  considered  in  the  weight  equations.  Using 
this  assumption,  the  weight  solutions  for  a two-element  adaptive  array 
were  derived  by  reducing  the  order  of  the  equations  describing  the  sys- 
tem. An  analytical  method  for  determinign  the  weight  values  with  a 
continuous  spatial  distribution  of  interference  is  also  presented. 

From  the  weight  solutions,  we  have  observed  that  the  element 
noise  enables  the  jamming  signal  to  be  present  at  the  array  output 
when  the  array  is  fully  constrained.  When  the  array  is  over- 
constrained, the  interference  is  always  present  at  the  array  output. 

Also,  control  of  the  desired  signal  phase  may  be  lost  when  many 
desired  signals  are  present. 

With  a continuous  interference  sector  present,  we  have  seen 
that  the  SNR  degrades  as  this  sector  increases,  when  the  desired 
signal  is  separated  from  the  interference  by  more  than  7°.  when 
this  separation  becomes  small  (<7°),  we  find  that  the  SNR  improves 
as  the  jamming  sector  width  increases. 
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APPENDIX  I 

In  this  appendix  we  discuss  a method  for  determining  the 
coordinate  rotation  matrix  R given  in  Eq.  ( 57)  from  its  reduced 
form  R given  by  Eq.  (100). 

The  equations  describing  the  response  of  a two-element  adaptive 
array  are  given  in  Eq.  (57).  It  was  noted  by  Compton  |21|  that  when 
: exhibits  the  type  of  symmetry  shown  in  Eq.  (58),  the  four  coupled 
equations  of  Eq.  (57)  may  be  reduced  to  a pair  of  complex  equations 
and  the  eigenvalues  of  may  be  determined  from  the  reduced  form  of 
■>  . It  will  be  shown  in  this  appendix  that  the  eigenvectors  and 
hence  the  coordinate  rotation  matrix  may  also  be  found  from  the 
reduced  equations. 

It  was  noted  in  section  III  that  by  allowing 
( A- 1 ) A = *n  + j*12 
(A-2)  B = *13  + 

(A- 3)  C = 4*  3 j + j4l32 
(A-4)  D = *33  + j«3l+ 
a new  matrix  Q can  be  formed  given  by 


87 


(A-5)  Q = 


A B 


C D 


Then  by  finding  the  eigenvalues  of  Q,  the  eigenvalues  of  # may  be 
determined  |2l|.  Specifically,  if  the  two  eigenvalues  of  Q are 


then  the  eigenvalues  of  * are 


(A-7)  = Ll.  L*.  L , It 


Recall  that  it  was  also  shown  in  section  II  that  the  equation 
of  a two-element  array  may  be  reduced  to  the  form 


(A.8)  - (v**  U =2k 

l W3-jW4  J V^3]+J^32  ^33+^34  AW3-^W4  / 


Sl-J$2 

S3‘JS4 


In  order  to  solve  Eq.  (A-8)  we  make  a rotation  of  coordinates.  Let 


(A-9)  W = 


wi  - jw2\  / r12V* 


W3  - JW 


?2i  r22/ty 


™ m 


(A-6 ) Aq  • ± Lj,  ^2  ± _ L2 


-if 


Substituting  Eq.  (A-9)  into  Eq.  (A-8)  and  multiplying  the 
result  on  the  left  by  R-1  gives 


(A-10) 


d_i 

dt 


+ 2k 


v 2 / 


(}n  f 12)  (*  “V1 


% 'V/ 

yr21  r22/ 


B r 


\ 


11  ‘12 


V 


a,  K 

C D/\r 


2 1 r22 


■y 


(*  \ 
1 yi 

r\j 

= 2k 

. r\  1 

\y  j 

1 

CM 

tH 

CM 

>c_ 

where 


(A-ll) 


fs  1 ' J 


s -J  s4 


If  the  rotation  matrix  R is  chosen  so  the  bracketed  matrix 
product  is  diagonal. 


(A-12) 


[ ? ? 

11  12 


V1 

A 

B ' 

( r 

11 

r ^ 
12 

0 ^ 

1 

■v 

*\t 

/ 

Vc 

r22  y 

w 

then  Eq.  (A-8)  are  uncoupled  and  the  solution  is  easily  obtained. 

The  values  1^  and  L2  are  the  eigenvalues  of  Q.  The  unitary  rotation 
matrix  R is  composed  of  the  eigenvectors  E2  and  E2  corresponding  to 
the  eigenvalues  Lx  and  L2, 
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r*  l/>  ? 10 


AD-A036  581 
UNCLASSIFIED 

2 or  2 

AD 

A03658I 


OHIO  STATE  UNI V COLUMBUS  ELECTROSCIENCE  LAB  F/G  9/5 

ADAPTIVE  ARRAY  PERFORMANCE  WITH  SPATIALLY  DISPERSED  INTERFERED— ETC  (U) 
JUN  74  D M D1  CARLO  N00019-74-C-0141 

ESL-3832-1  NL 


Since  R is  unitary 


(A- 14)  fT1  = Rt  . 


We  can  then  rewrite  Eq.  (A-12)  as 
(A— 1 5 ) 

Xii-Jyn  x2i“Jy2iV A 


'v 


B \(  xu+jyn  *12+Jy  \ 


'Vx 


C12"jy!2  X22-jy22AC  D 


(A-16)  r. . = 

1J  JT 


From  Eq.  (47)  we  have 


( A- 17) 


All  r2i  r 3 1 r41  \4l]  *12 


RT*R  = 1] 


r 12  ril  F 32  F 42 


l + j *■  0 

1 2 


VX2  1+Jy21  X22+Jy22y 

where  the  elements  of  R are  now  defined  as 


0 * 


F 1 3 f23  r 3 3 r43 


\^14  r24  F 34  F 44/  *41 


13  *14  ^ f 

21  *22  4*2  3 *24 

31  4*  32  4>33  <f>34 

42  4*  43  4>44 


rH  ri2  f13 


r r r 

21  1 22  *23 


r 31  r32  F 3 3 


r r r 

41  *42  J43 


! s 0 0 

-v2  t, 00 

! 0 0 *3 

\ o o 

V d 


(See  | 21]  pp.  90-91.) 

From  Eq.  (A-17)  we  notice 


(A-18)  "e»l(L|).i1.jL  ji  j)ri,njrji 


(A-19)  In»g[L,]-s4  j,  j,  ru*ij  rj2- 


The  following  relations  are  evident  from  the  symmetry  shown 


in  Eq.  (58): 


(A-20a)  <j>jj  = <P22 


(A-20b)  *l2  = <t>2l 


(A-20c)  3-  ^2l4 
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(A-20d) 

*14  = "*23 

(A-20e) 

*31  = *42 

(A-20f) 

*32  = "*41 

(A-20g) 

*33  " *44 

(A-20h) 

*34  = "*43 

Now  if  we  expand  Eqs. 
and  (A-19)  and  make  a 
ships  of  Eqs.  (A-20), 

(A-15)  and  (A-17)  into  the  form  of  Eqs.  (A-18) 
term-by-term  comparison  imposing  the  relation- 
we  find 

(A-21a) 

ril  = X11 

( A-2 lb ) 

r21  = 1 1 

(A-21c) 

r 31  = X2  1 

(A-21d) 

F41  = ”^21 

(A-22a) 

ri2  = ^11 

(A-22b)  r 


(A-22c)  r 


(A-22d)  r 


(A-23a)  r 


(A-23b)  r 


(A-23c)  r 


(A-24b) 


(A-24c)  r 


(A-24d)  r 


In  other  words 


I 


where  and  are  the  term  of  R given  in  Eq.  (A-16).  Then  by 
finding  the  complex  rotation  matrix  R,  we  have  also  determined  the 
real  rotation  matrix  R.  This  is  also  to  say  that  by  finding  the 
eigenvectors  of  Q,  we  have  found  the  eigenvectors  of  $ . 
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1 nJMENSlON  AL PH  < 5 ) 

2 roMPLt*  AlPH 

3 r********PATTrRN  ROUTINE********** 


4 

nlPENSlOiM  A ( 1 2 ) « T HET  A ( 12 ) , SIGU2).  SIPJ12) 

5 

COMPLEX  TEN.  XYZ 

6 

COMPLEX  LESO,  ftjvc 

7 

COMPLEX  XSJ,  X c , HjC  ,W2C 

8 

COMPLEX  SlF,X,T,U,SRP,XX,TT,IIU 

9 

COMPLEX  YY.ZZ»SS,  AAN 

10 

OATA  PI.  1PI/3. 14159265,  6.2*31853/ 

11 

WRITE.  (8,61 > 

12 

61 

POKMAK ’TOTAL  mUMBER  OF  SIGNALS*) 

13 

RfAOCA,-)  N 

14 

WRITE ( 8 *77) 

15 

77 

FORMAT ( * TOTAL  NUMBER  OF  OESTREO  SIGNALS*) 

16 

Rf  AO ( 6 » - ) M 

17 

WRITE  (f *fc2 ) 

16 

62 

r0RMAT(  ‘SIGNAL  LOCATIONS*) 

19 

no  2 J=1  , N 

20 

RE  AO ( 8 , « ) THTTA(J) 

21 

SIG( J)=l. 

22 

OUR=N«w 

23 

T Y=  THE ta(J)/( 57, 2957758) 

24 

M J)=PI*SIMTY) 

25 

2 

CONTINUE 

26 

WRITE (8,63) 

27 

63 

C ORMAT ( » NOISF  AMPLITUDE*) 

28 

PEA0<8,->  AN 

29 

nt  SO=CMPLX ( 0. ,0,  ) 

30 

AjMC=CvPLX(0.,n.) 

31 

onTsO  . 

32 

nu  21  1=1, N 

33 

no  21  J=1 , N 

34 

?1 

RNT=RNT+(1.-COS(A(I)-a(J) ) ) 

35 

ysCMPLV (U, , 0.  ) 

36 

y s=  x 

37 

no  5 K=l , N 

36 

WKI TE ( 6 . • ) THFTA|«) 

39 

XX  = CMPl  X ( 0, »*A  <K ) ) 

40 

*SJ=CMPLX(0. « A (K)  ) 

41 

XS=XS*SIG(R )*STG(R)*CEXP(XX) 

42 

5 

c UN  TIN' IE 

43 

44 

4 5 

46 

47 
40 

49 

50 

51 

52 

53 

54 

55 

56 

57 
56 

59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 
83 
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ru  7 k-i,*, 

*x=c“iplx ({.. , *M*n 

7 V = X + SI<*-<K)*CFXP(XX) 

Y=n. 

HO  8 

v = Y4Slr.(K)*sir,iK) 

7 = 0. 

rONTIMi'K 
P roNTT'M*  *■! 

Y=  Y+2.  ♦ AAi*A*. 
n0  9 K = l»* 
n 7 = 7 *S  1 ( K ) 
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6 = 0 . 
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niMENSlON  I RUF (100) 

HAT  A Pi, TPI, OPR /3. 14159265, 6. 2831 653. 57. >957795 
r aLL  PLOTS( IBUF»10C,3) 

TALL  FACTOR ( ,5 ) 

CMXsO* 

I 1 = 2 
I 0 = LI-i 

UNITE (1,777)  IPFO.  ISS 
777  roRWAT ( • 1PFQ=»,I>,  • ISSs'.t?//) 
nO  29  Jsl, 361,1 
t iELQsfO.  ,0.  ) 
no  30  IJsl.LI.l 
TU=IJ-1 
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AK«  ( I J>  = ' I * C 1 . / I . * *ros  ( J/l-MP  1 

ft. 

"'.A 

akT  { IJ'SCH'I  V ( n.  ./  4k(  1.J)  ) 
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P 1(  LfsP  If.  LO  + Al  PH(  I J ) *Cf  XR  ( in*n<T  ( T J)  > 

9b 

?e 
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n 
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97 
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it 

99 

24 

COMTIM-F 
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in  J 
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102 

f i\if' 

i 

i j 

103 
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‘i 

10“ 

r i vt<v,bi  o<*  cf  ( . p,i  ) 

1 Ob 
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|j 

10. 
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i 

1"7 

r’K  110  T=1.o 

If  6 

BflOrf'  . ♦ < 5-  J > /4  . 

h 

104 

r ,.L  L tn.nT  (0/  D.o.  , 3) 

1 

110 

ro  11'  Jrl.ifi 

111 

AAjG=<  J-l >»PJ/9P. 

to 

112 

y XzR  AD*COjs  ( ANG) 

113 

yy=RAO*Sll'!(  ONG  ) 

4 

114 

110 

tall  plot (xy.yy.2 ) 

115 

no  112  i=i.t. 

tt 

116 

DhGsI  I.l ) .1  I/O. 

k 

117 

AlvpS=ANf;*lJI 

lie 

A t'.  G F z A f'  G 

119 

TF<  I .(•  ,2*< 1/2)  > GO  TO  111 

i 

120 

A >Y  6 S z A ■ 0 

121 

AlviGFzA  'G  + PT 

1?? 
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roMTTMif 

I 

1 23 

yx=5. *f05 ( A'  09 ) 

» 

124 
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125 
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to 
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129 
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to 

130 
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131 
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r 

132 
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I 

133 
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u 

134 

YY=RAO*SlN<  AA<G1 

135 

XX=RAU*COS<  A,\G> 

* *! 
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136 
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• 

1 

137 
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3 2n 
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n 

139 
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1 

140 
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141 
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rtvn$ 

1 
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